Recombination of three ultracold fermionic atoms 
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Three-body recombination of identical, spin-polarized fermionic atoms in the ultracold limit is 
investigated. The mechanisms for recombination are described in terms of the "scattering volume" 
V p in the framework of the adiabatic hyperspherical representation. We have calculated numerically 
the recombination rate Kz as a function of V p and have found that Kz scales as |Vp| 8//3 for small 
\Vp\. A comparison with experimental data is also presented. 
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Recently, the quantum degenerate regime was attained 
in ultracold gases of fermionic atoms such as 40 K and 
6 Li [Q, ||. Part of the motivation for these experiments is 
to observe a pairing of fermions, leading to a superfluid 
state. One important factor limiting the achievable den- 
sity in these degenerate Fermi gases (DFG's) of trapped 
atoms is the loss of atoms through three-body recom- 
bination. Such losses occur when three atoms scatter 
to form a molecular bound state and a third atom - 
K + K + K — > K2 + K, for instance. The kinetic energy of 
the final state particles causes them to escape from the 
trapping potential. 

While ultracold three-body recombination of identical, 
spin-polarized bosons has been theoretically studied be- 
cause of its importance for Bose-Einstein condensates, 
recombination of identical fermions has not yet been con- 
sidered. For bosons, Fedichev et al. B predicted that the 
recombination rate K3 grows with the two-body s-wave 
scattering length a s , namely K 3 oc a 4 , for a s > 0. This 
scaling was later confirmed by Nielsen and Macek |B| who 
also pointed out that it should hold for negative a s . The 
a 4 scaling law for both signs of a s was indeed obtained 
by Esry et al. ||, Bedaque et al. |Q, and Braaten and 
Hammer j^]. 

In the case of fermions, however, the Pauli exclusion 
principle prohibits s-wave scattering of atoms in identi- 
cal spin states, thus leaving only p-wave collisions. The 
relevant low-energy scattering parameter in this case is 
the two-body p-wave "scattering volume" defined as 
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where 5 p (k) is the p-wave scattering phase shift and k is 
the wave number. The scattering volume V p is related to 
the p-wave scattering length a p (see, for instance, Ref. H) 



by V p 



We choose V p , rather than dp , ctS the pa- 



rameter to characterize the three-body recombination of 
fermions since an artificial nonanalyticity is introduced 
into a p when taking the cube root of the quantity in the 
right-hand side of Eq. (Q). 

Even though recombination of identical fermions is 
suppressed at ultracold temperatures by the Pauli princi- 



ple, it does not vanish. In fact, it has been shown that the 
rate is proportional to E 2 at low collision energies Jl2| . 
While this rate remains negligible under typical experi- 
mental conditions, it can become substantial near a Fesh- 
bach resonance. The E 2 threshold law no longer applies, 
and the recombination rate tends to the limit imposed by 
unitarity — often comparable to or larger than the rates 
for boson systems. Feshbach resonances are, of course, 
extremely useful tools for the experimentalist, so under- 
standing the behavior near such a resonance is crucial. 
So far, such resonances have been observed, for exam- 
ple, in systems of 40 K Jl(J and 6 Li This resonant 
tuning was also observed for other alkali species in BEC 
experiments |Tx| . 

This Letter treats the three-body recombination of 
identical, spin-polarized fermions in the ultracold limit. 
We examine the recombination rate K3 as a function 
of the scattering volume V p by numerically solving the 
three-body scattering problem. First of all, since the 
atoms are spin polarized and thus in a completely sym- 
metric spin state, the spatial wave function must be com- 
pletely antisymmetric in order to satisfy the Pauli prin- 
ciple. This fact, combined with the generalization of 
Wigner's threshold law to K3, shows that the J n = 1 + 
symmetry dominates at threshold, where J is the total 
orbital angular momentum and n is the parity of the sys- 
tem. It follows that the recombination rate depends on 
the collision energy E as E 2 near threshold p2[ . (The 
same analysis applied to the boson case yields the fa- 
miliar result that the recombination rate is constant at 
threshold due to the + symmetry.) Therefore, we con- 
sider only the 1 + case. Dimensional analysis, together 
with this E 2 law, suggests a |V^| 8 / 3 scaling of the recom- 
bination rate K 3 as opposed to the a 4 scaling for bosons. 

The interaction potential used is a sum of triplet two- 
body potentials, i.e. V = v(ri 2 ) + v(r 2 3) + v(r 3 i). This 
choice is appropriate for fully spin-polarized atoms that 
collide in a quartet electronic state. For simplicity, we 
model the two-body potential as either 
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or 
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The former potential has proven convenient in recombi- 
nation calculations while the latter has a more physical 
van der Waal's tail with Cq = r^D. Ideally, the results 
will not depend on the particular potential used in the ul- 
tracold limit. The parameter ro controls the range of the 
potential. The coefficient D, representing the potential 
depth, is treated as an adjustable parameter that per- 
mits us to control the scattering volume V p , thus mim- 
icking the tuning ability of an external magnetic field. 
This coefficient also economizes the calculations substan- 
tially since it allows us to reduce the number of two-body 
bound states without sacrificing the ultracold physics. 

The two-body p-wave scattering volume V p behaves 
much like the s-wave scattering length as a function of the 
parameter D, displaying a tangent-like structure. As D 
becomes more negative and the potential becomes more 
attractive, the scattering volume passes through a pole 
and changes sign each time the potential becomes deep 
enough to support one additional p-wave bound state. 
For simplicity, we consider the parameter range for which 
there exists only one two-body bound state. 

The details of the theoretical methods we employ are 
more completely discussed in Ref . Q , so only a brief out- 
line is given here. We use the adiabatic hyperspherical 
representation with modified Smith- Whitten hyperspher- 
ical coordinates. Simply put, hyperspherical coordinates 
transform the six relative Cartesian coordinates in the 
center of mass frame to a set with a single length co- 
ordinate, the hyperradius R, and five hyperangles. The 
hyperradius can thus be thought of as characterizing the 
overall size of the three-body system. These coordinates 
also allow us to easily impose the correct permutation 
symmetry on the wave functions. Solution of the adia- 
batic equation yields adiabatic hyperspherical potential 
curves and channel functions. The coupled hyperradial 
equations are then solved using an i?-matrix propagation 
method. 

The adiabatic hyperspherical representation reduces 
the collision of the three atoms to dynamics on a set 
of coupled hyperradial adiabatic potentials. These po- 
tentials bear a strong resemblance to standard molecular 
Born-Oppenheimer potentials (see Fig. and can be in- 
terpreted in much the same way. The lowest potential 
curve in the figure, for instance, correlates to a bound 
molecule and a free atom far away. Because the atoms 
are identical fermions, the molecule has unit angular mo- 
mentum and the potential has a centrifugal barrier. Fur- 
ther, because we consider only the f + symmetry due to 
its dominance at threshold, the free atom must also be in 
a p-wave relative to the center of mass of the molecule. 
All of the other curves in Fig. fil correlate to three free 
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FIG. 1: The lowest 12 adiabatic hypersperical potential 
curves for three spin-polarized fermions with V p = —1.25 x 
10 5 a.u. 3 using the sech 2 potential with ro=15 a.u. 

atoms. In fact, there are an infinite number of potential 
curves associated with three-body continuum channels 
and approach the three-body breakup threshold {7 = 
asymptotically. 

Three-body recombination can be seen as a transition 
from one three-body continuum channel to the recombi- 
nation channel, driven by nonadiabatic coupling. Qual- 
itatively, the adiabatic potentials shown in Fig. [l] dis- 
play analogous behavior, as functions of V p , as the adia- 
batic potentials for bosons show as functions of a s H [Tt]] 
For both systems, the entrance channel goes from be- 
ing strongly repulsive for positive V p and a s to having 
an attractive well behind a potential barrier for nega- 
tive values. The entrance channel for the fermions is 
generally more repulsive, however, than the addition of 
a simple J=l centrifugal potential term to the boson 
curves would give due to permutation symmetry consid- 
erations p2| . The recombination channels are also very 
similar, although it should be noted that there is no Efi- 
mov effect [|l4| for fermions in the limit \V P \ — » 00. Efi- 
mov physics plays a key role in the interpretation of the 
ultracold recombination of bosons |^| . 

The primary difference between fermion and boson sys- 
tems lies in the nonadiabatic coupling. While it is still 
similar for negative V p and a s , for positive values the 
similarities end. Where the coupling strength for bosons 
shows a definite peak whose position increases linearly 
in a s , the coupling strength for fermions remains peaked 

at small R with a slowly decaying shoulder whose ex- 

1/3 

tent grows in proportion to V p . For negative values of 
V p , we thus expect that the fermion recombination rate 
can show resonant enhancement due to three-body shape 
resonances just as for bosons, but we do not expect an 
infinite series of such resonances since this was a conse- 
quence of the Efimov effect. We have not, however, seen 
evidence of such resonances as V p — > — 00 in our calcula- 
tions. 

To calculate the three-body recombination rate, we 
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solve the hyperradial equations |T]| for potentials like 
those shown in Fig. |l|. The numerical hyperradial wave 
functions are matched to the appropriate Bessel func- 
tions based upon the known form of the long-range poten- 
tials |l3j. The iS-matrix is obtained and the generalized 
cross section for three-body recombination is calculated: 
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Here, k = ^2/j J E/h 2 is the hyperspherical wave number 
in the incident three-body continuum channel, and the 
indices A and / label initial three-body continuum and 
final recombination channels, respectively. As it turns 
out, the numerical prefactor in Eq. ([|) — which is de- 
termined by permutation symmetry — is the same as for 
three identical bosons j(| (recall that we are considering 
the = 1 + symmetry). The event rate constant per 
atomic triad is then defined simply as K 3 = ^cr. This 
quantity is related to the atom-loss rate constant L 3 by 
L 3 = 3K 3 /6 @. 

Because the rate depends strongly on the collision en- 
ergy in the ultracold regime, it must be thermally aver- 
aged in order to compare with experimental data. Fol- 
lowing Ref . Jill , we have derived the thermally averaged 
recombination rate constant to be 
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We show in Fig. || (K 3 ) as a function of V p for a temper- 
ature of 2 fiK. Specifically, Fig. | shows (K 3 /rffi/ 8 versus 
V p /tq. We plot K^ 8 to best reveal whether the predicted 
|VJ,| 8 / 3 scaling actually holds for the calculated rate; the 
factors of ro were included based on general arguments 
of length scale invariance of the Schrodinger equation. 
Because the two-body potentials in Eqs. (||)-(|3|) depend 
only on the combination ry/ro, ro can be defined to be 
the new length scale for both the two-body and three- 
body systems. This new length scale implies multiplying 
energies and temperatures by Tq, lengths by r^ 1 , and the 
rate by Tq . The energy parameter D in our model — 
which controls V p — is independent of the length scale, 
and so its scaled counterpart, r^D, can be used to control 
V p /tq. Having calculated the rate as a function of V p and 
energy for a fixed ro, we can thus obtain the rate for any 
other value of ro by means of the above scaling. The key 
is that at a fixed laboratory collision energy (or temper- 
ature), changing ro effectively changes the range of the 
two-body potential, and the new rate corresponds to the 
rate at the scaled energy r^E. One has to be careful, 
of course, to track the simultaneous change in the scaled 




FIG. 2: The recombination rate K3 as a function of V p for 
potentials from Eq. (Q) (no symbols) and Eq. (^|) (diamonds) . 
The solid lines corresponds to the thermally averaged rate 
at T=2 [iK; and the dashed lines, to the energy-dependent 
rate at E=2 /iK. All curves were calculated with ro = 15 a.u. 
The inset focuses on the region with small \V P \ showing a 
Stiickelberg oscillation. 



The solid lines in the figure show the thermally av- 
eraged rates (K 3 )(T), while the dashed lines show the 
energy-dependent rates K 3 (E). The symbols for the 
sech 2 potential, Eq. (||), are omitted for clarity since the 
rate was obtained on a relatively dense set of V p values. 
The diamonds denote the rates resulting from the second 
potential, Eq. (^). The rates from the two potentials — 
which are qualitatively very different — show reasonable 
quantitative agreement, especially for the more experi- 
mentally relevant thermally averaged rate. The poorest 
agreement occurs at small \V P \ where the rates are also 
small. If, however, the rate curve for either potential with 
a different value of ro and the same lab energy E were 
plotted here, it would look systematically different, al- 
though qualitatively similar. This apparent dependence 
on the two-body potential would seem to argue against 
a "universal" curve as was found for bosons J|, |6[ 0, ||. 
Based on the above length scaling arguments, though, it 
is clear that the scaled rate curves — at the same scaled 
energy — are what should be compared between differ- 
ent potentials and would thus come closest to a universal 
curve. 

Figure H shows that the |V^,| 8 / 3 scaling does indeed hold 
quite well for less than roughly 2500 for K 3 (E), 

but for (K 3 )(T) it only holds over the region | ] / r § be- 
low approximately 800. Since this scaling is based on the 
threshold behavior of the rate, it is not surprising that 
it breaks down sooner for the thermally averaged rate 
given that the averaging procedure includes rates from 
higher energies. The fact that the scaling breaks down 
for K 3 (E) is an indication that the fixed collision energy 
is no longer in the threshold regime. 

For large positive scattering volumes V p , the collision 
energy eventually becomes large compared to the dimer 
binding energy, and K 3 increases more slowly than the 
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|T<p| 8 / 3 scaling law. Recalling that for negative V p there is 
a barrier in the initial adiabatic hyperspherical potential, 
the breakdown in the scaling law comes when the barrier 
sinks below the collision energy at large negative V p . The 
fact that the recombination rate increases more rapidly 
than the scaling law with increasing \V P \ simply reflects 
the more ready tunneling through the barrier. The rate 
continues to increase until it reaches the limit imposed 
by unitarity, namely (in SI units) 
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This value is obtained from Eq. by assuming that 
the recombination predominantly comes from the lowest 
continuum channel with unit probability. We expect that 
this limit will be reached when \V P \ becomes large enough 
that the barrier in the three-body entrance channel is 
lower than the collision energy. For the conditions of 
Fig. ||, this value of V p is well off the scale of the plot. 

Interestingly, the inset of Fig. || shows a small peak in 
the rate located at about Vp/r^ — —60, i.e. at B=200.8 
G. This peak is a constructive interference between two 
indistinguishable pathways — a Stiickelberg oscillation. 
This peak disappears and reappears as a function of col- 
lision energy since the phase difference between the two 
pathways depends on the energy. 

In Fig. ||, we compare our calculations with experi- 
mental measurements of the recombination rate of three 
40 K atoms all in the 1 9/2 , —7/2) spin state near ap-wave 
two-body Feshbach resonance. The theoretical rates are 
calculated at T = 2.5 /iK, while the experimental tem- 
perature varies between 2 and 3 ^iK. The figure shows 
the rates as functions of the magnetic field strength B. 
The Feshbach resonance lies at Bq — 198.5 G; V p is posi- 
tive for B below the resonance and negative for B above. 
Treating ro as a fitting parameter, we obtain the solid line 
shown in the figure. The agreement on the rise of the res- 
onance (V p >0) is better than on the fall (V p <0). The 
disagreement away from the resonance is likely due to 
the breakdown of the experimental assumption of purely 
three-body loss in determining the rate coefficient. Near 
the resonance, however, our theoretical rates between one 
and two orders of magnitude lower. In fact, for all values 
of ro, the theoretical values lie below the experimental 
ones. This discrepancy is not surprising given that we 
included only a single partial wave (1 + ) in our calcula- 
tion. While this symmetry controls the limiting behav- 
ior at threshold, it is not necessarily dominant near the 
resonant peak. In fact, limited explorations have demon- 
strated that other J 77 can be comparable to 1 + at the res- 
onance peak. The contributions from these higher partial 
waves accumulate quickly since thermal averaging, from 
Eq. (||), gives an extra factor of roughly (J + 2)!. 

In summary, we have investigated ultracold three-body 
recombination of identical spin-polarized fermions and 
have numerically calculated the recombination rate as a 
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FIG. 3: The theoretical rate (solid line) and the rate mea- 
sured by C. Regal et al. (filled circles) are plotted as functions 
of the magnetic field strength B. 



function of the scattering volume. This is the first such 
study that we are aware of, and much remains to be un- 
derstood about these fermion systems. Moreover, with 
some slight modifications to our numerical method, it is 
possible to treat systems with only two identical parti- 
cles. Three-body recombination of such systems plays an 
important role in experiments on mixed-spin state Fermi 
gases and mixed Bose-Fermi gases. 
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